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Abstract 


We apply a primal-dual simplex algorithm for solving the biobjective min- 
imum cost-time network flow problem such that the total shipping cost and 
the total shipping fixed time are considered as the first and second objective 
functions, respectively. To convert the proposed model into a single-objective 
parametric one, the weighted sum scalarization technique is commonly used. 
This problem is a mixed-integer programming, which the decision variables 
are directly dependent together. Generally, the previous works have consid- 
ered the linear biobjective problem with the traditional network flow con- 
straints, while in this paper, corresponding to each flow variable, a binary 
variable is defined. These zero-one variables are utilized to describe a fixed 
shipping time for positive flows. The proposed method is successful in finding 
all supported efficient solutions of a real numerical example. 
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1 Introduction 


In multi-objective network flow optimization, several objective functions, 
which generally are time, cost, risk, environmental concerns, and so on, have 
to be optimized such that the flow conservation and capacity constraints are 
satisfied. In actual problems, since the objectives are frequently conflicting, 
finding a feasible solution optimizing all objectives simultaneously is difficult 
or in most cases is impossible. Therefore, the aim of multi-objective opti- 
mization is to compute nondominated points rather than optimal values. A 
perfect review of the multi-objective cost flow problem is done in [14, 25]. 

The minimum cost-time network flow problem (MCTNFP) is a well- 
known and valuable topic in the scientific researches, which can be classified 
in a subset of the fixed charge network flow problems (FCNFPs). The first 
objective function considers the total shipping cost of flow, and the second 
objective function measures the total fixed shipping time of flow. Correspond- 
ing to any flow variables, there is a zero-one decision variable that is utilized 
to compute the total fixed shipping time on each directed arc. The weighted 
sum scalarization technique is used to convert the reformulated model into a 
parametric single objective optimization. Eusébio, Figueria, and Ehrgott [11] 
presented a modified primal-dual simplex algorithm based on the work done 
in [8] to find nondominated extreme points of the linear biobjective network 
flow problems. The main idea in both algorithms is that the scalarization 
parameter \ € (0,1) is broken into subintervals such that each nondominated 
extreme point is associated with one and only one set in the partition. Ac- 
cording to these papers, the primal-dual approach is considered for computing 
all efficient solutions of MCTNFP (see also [22, 4]). 

FCNFP is a classical NP-hard combinatorial problem, (see [13]). This 
problem consists of two costs corresponding to each incident arc in the net- 
work: A fixed cost for the use of the arc and a variable cost proportional to 
the number of units sending along the arc (see [7]). The end of the FCNFP 
is to choose some arcs and assign feasible flows to them in order to transfer 
commodities from an origin to a destination at a minimal total cost. Some 
important traditional problems such as transportation problem, facility loca- 
tion problem, and network design problem can be formulated as an FCNFP. 
Many exact and heuristic techniques have been applied for solving the FC- 
NFP. Analytic procedures commonly utilize branch and bound method to 
search an exact solution of the FCNFP (see [21]). Heuristic and evolutionary 
techniques have been developed to find a near-optimal solution to the under 
regarded problem (see [19, 20]). 

Raith and Ehrgott [24] presented an algorithm to compute a complete set 
of efficient solutions for the biobjective integer minimum cost flow (BOIMCF) 
problem based on a two-phase method (see [23]). Using the mean of weighted 
sum technique, Abo-Sinna and Rizk-Allah [2] introduced a new gradient- 
based neural network approach for solving biobjective optimization, which 
is stable in the sense of Lyapunov. According to Benson’s scalarization [6] 
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and min-max methods, Mohammadi, PourKarimi, and Pedram [17] stud- 
ied the scientific workflow scheduling problem in a multi-cloud environment 
with cost and makespan minimization objectives. Eusébio and Figueria [10] 
solved a sequence of e-constraint problem to find all integer nondominated 
extreme points of small or medium size of BOIMCF problem by a branch- 
and-bound method. Ehrgott, Shao, and Schébel [9] implemented a method 
for approximating the nondominated set of a multi-objective nonlinear con- 
vex programming problem. Abareshi and Zaferanieh [1] introduced a bilevel 
model to solve the capacitated p-median facility location problem with the 
most likely allocation solution. Applying Lagrangian dual theory, the pro- 
posed bilevel problem is reduced to a new one-level nonlinear mixed-integer 
problem whose solution is obtained by comparing two mixed-integer linear 
problems. 


In this study, we introduce a strategy for formulating and solving biob- 
jective minimum cost-time network flow problems (BOMCTNFPs), which 
in addition to the flow variables, some corresponding zero-one decision vari- 
ables are used in their construction. Keshavarz and Toloo [16] considered 
the weighted sum scalarization approach to convert their proposed model to 
a corresponding parametric mixed-integer minimum cost flow problem with 
single objective function and finally reported the supported efficient solutions 
for only two values of the scalarization parameter \ = 0.05, 0.5. The main 
difference of our work with [16] is that, by utilizing the primal-dual technique, 
all of the nondominated extreme points and the supported efficient solutions 
of the proposed problem are computed. 

The idea of the mathematical formulation of BOMCTNFP is taken from 
Hochbaum and Segev [15]. They offered an augmented formulation of the 
fixed charge problem and proved that an optimal solution of the augmented 
problem is an optimal solution of the fixed charge problem. The numerical 
procedure that they have been implemented is based on Lagrangian relax- 
ation, which behaves well empirically. 


2 Mathematical formulation of the BOMCTNFP 


Ehrgott and Puerto [8] suggested the primal-dual algorithm for solving a gen- 
eral form of the multi-objective programming problem. This algorithm solves 
the problem by constructing a partition of the set (0,1) such that each subin- 
terval of this partition is attributed to only one efficient solution. A modified 
primal-dual approach for solving linear biobjective minimum cost network 
flow problem was developed in [11], which uses reduced cost information to 
avoid redundancy. Here, we apply a similar idea to find nondominated points 
for solving a mixed BOMCTNFP. 

Suppose that G = (V,.A) is a directed and connected network in which V 
is a finite set of nods with |V| = m, and A is a collection of directed pairs of 
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elements of V called arcs, with |A| =n. Corresponding to any arc (i,j) € A, 
a triple label as (c¢;;,t;;, uj;) is assigned, where the first, second, and third 
components are shipping cost, shipping fixed time, and capacity, respectively. 
Associated with every node i € VY, there is a value b;, which shows the supply 
(if b; > 0) or demand (if b; < 0). If b; = 0 for some 7 € V, then the node i is 
called a transshipment node. The interested reader can consult [5, 3]. 

Consider the following BOMCTNFP with mixed continuous-binary vari- 
ables: 


mincr = ) Cig Liz, 


(i,j)EA 
minty = S- lig Yig, 
(A,j)EA 
Bat SY) ey — D> tin = be for all k € V, 
(k,j)EA (i,k)EA 
0< Lij < Uij for all (i, 7) E A, (1) 
ro 1 Liz > 0, 
va 0 Lig = 0, 
Yj © {0, 1} for all (i, 7) EA, 


where the binary variable y;; is an auxiliary variable to consider the fixed 
time t;; for the positive flow x;; on arc (i,j). If we assume that the values of 
the flow variables are real, then the model (1) is a biobjective mixed-integer 
problem, and if the value of the flow variables is integer, then this model 
is a biobjective integer problem. We refer to the first case with minimum 
cost-time continuous flow (MCTCF) problem, and to the second one with 
minimum cost-time integer flow (MCTIF) problem. 

The efficient solution is an important concept in the problem of multi- 
objective optimization. In fact, an efficient solution to a multi-objective 
problem is the response, which cannot improve some of the objectives with- 
out worsening other objectives. Assume that X is a set of feasible solutions 
to problem (1), providing the constraint y;; € {0,1}, (i,7) € A is replaced 
with the continuous constraint y;; € [0,1], (i,j) € A. The set X is called 
the decision space. We denote the feasible set in the objective space as 
Z(x,y) = {(cx, ty)|(w, y) € X}. 


Definition 1. [12] Let (x,y), (a’,y/) © X. If cx < ca’,ty < ty’ and 
(ca, ty) # (ca’, ty’), then it is called (x,y) dominates (a’, y’) in the decision 
space, and equivalently (cx,ty) dominates (cx’, ty’) in the objective space. 
This notion is denoted by Z(a,y) x Z(2’,y’). 


Definition 2. [18] The feasible solution (x,y) € X is an efficient so- 
lution, if there is no other feasible solution like (x’,y’) € X such that 
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Z(a',y') ~ Z(ax,y). If (x,y) is an efficient solution in the decision space 
X, then its image in the objective space is non-dominated, that is, Z(z, y) is 
named a nondominated point. The set of all efficient solution is denoted by 
Xf, and the set of all nondominated points is denoted by Zp. 


Definition 3. [12] The efficient solution (x*,y*) € Xz is a supported effi- 
cient solution, whenever it is an optimal solution to the following weighted 
objective problem: 


min{A;cx + Agty: (x,y) € Xf, (2) 
where (Ai, A2) € Ap and Ap = {(A1, A2) : Ay > 0, A2 > 0, Ay t+AQ = Lt, 


Remark 1. [8] If (2*,y*) is a supported efficient solution, then Z(x*, y*) is 
also called a supported nondominated point. We use X,~ and Zn to denote 
sets of supported efficient solutions and supported nondominated points, re- 
spectively. The efficient solution (2*, y*) is a nonsupported efficient solution, 
if there exist some nonpositive values of Ay or Ag such that (x*,y*) is an 
optimal solution of problem (2). 


3 An augmented formulation of the BOMCTNFP 


In this section, we look for the supported and nonsupported efficient solutions 
that are the main challenges in solving mixed BOMCTNEFPs. In addition, we 
show the relation between x;; and y;;, which is an additional challenge. We 
solve this problem with an alternative linear constraint, which is expressed in 
many articles in the field of fixed-charge problem (see [15, 7, 20]). First, we 
reformulate model (1) by replacing the third constraint with 2; < uijyi,;, 
as follows: 


mincr = S- CijLij, 
(i,7)€A 
minty = se lig Vig 
(i,9)€A 
s.t. S> try — D> win = dp forallkeV, (3) 
(k,j)E A= (i,k)EA 
Lig S Wiz Yay for all (1,7) € A, 
Liy > 0, Yay € {0, 1} for all (4,7) € A. 
Hochbaum and Segev [15] proved a theorem that states that an efficient so- 


lution to (3) is an efficient solution to (1). On the other hand, a fundamental 
theorem in multi-objective linear problem assures us that any efficient solu- 
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tion can be characterized as an optimal solution of the weighted sum single- 
objective problem (2). A proof for this fact can be found in [26]. Therefore, 
problem (3) can be reformulated as the following parametric problem, which 
provides a set of supported efficient solutions of the above problem: 


min ACL + (1 = A)ty = aN x Cig Lig + (1 = d) S- biG Vijs 


(4,j)EA (i,j)EA 
s.t S| oy So tees for all k € V, (4) 
(k,j)EA (i,kK)EA 
Lig < Wag Yay for all (t,7) € A, 
Liz =O, Yiy € {0, 1} for all (t,7) € A, 


for all A € (0,1). 


4 Primal-dual algorithm for solving mixed BOMCTNFP 


In the following, we use a primal-dual algorithm for the weighted sum single- 
objective FCNFP (4) and evaluate the efficiency of the algorithm. Since 
the primal-dual algorithm is based on the duality theory of problem (4), the 
relaxation relation y;; € [0,1] is used instead of y,;; € {0,1}. Hence, the dual 
problem of (4) is cast as follows: 


max So bate = S- Oi9; 


kev (i,j)EA 

st. 1 — 1} — pag < cij(A) for all (i,7)€ A, (5) 
[ij tig — O17 < taj (A) for all (i,7) € A, 
jig S 0, tg 0 for all (2,7) € A, 


where cj (A) = Aciz, tij(A) = (1— A)tiz, for some A € (0,1), = (™,.-., 7m) 
is the vector of dual variables corresponding to conservation constraint 
VEjer thi — Uiweatix = On, W = (Miz)G,je4 18 the vector of dual 
variables corresponding to the constraint xj; < ujjyij, and 6 = (61,5) G,j)ea is 


the the vector of dual variables associated with the constraint y;; < 1, that 
is, 


» Uj — Le Lig = dy : mt, forallk ey, 
(k,j)EA (i,K)EA 
Lig S Wag yig : Hi; for all (i,7) € A, (6) 


Yi < 13 bij for all (i, 7) € A. 
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Since the right side of the constraints in (5) depends on 4, the feasible 
solutions of this problem inevitably depends on \. Suppose that (27;, y;;) and 
(1 (A)*, u(A)*, 0(A)*) are the feasible solutions of the primal problem (4) and 
the dual problem (5), respectively. By using the complementary slackness 
conditions, (2*, y*) and (1*, w*,6*)(A) are optimized for their corresponding 
problems if and only if 


i) ( ai +t + ut 4 cig(A)) 2% = 0 for all i,j € V & for all (i,j) € A, 


ti) (- Mig tig + 65; — tiy(\)) vi =0 for all (,7) € A, 
iit) (- Ti, + Uiz¥iz ) Mig = 0 for all (i,7) € A, 
iv) (1 = uly) 5%, =0 for all (i,j) € A. 


If0 < aj, < wiz, then pF, = 0 and —77 +77 +ci;(A) = 0. Also, if ejj(A) — 77 + 
m; > 0, then 2, = 0 and y7, = 0. If cij(A) — 77 + 77 < 0, then consequently 
Hi; > 0 and xj; = yj;uij. Two cases may be occurred: 


a) If t;;(A) > ux,uij, then 6%, = 0; hence y¥, = 0 and a*, = 0. 
Hy ag 74d aw a A) 


b) If tij(A) < pijpuiz, then 67; > 0; hence y7; = 1 and x7; = wij. 


Let (7, ,5)(A) be a feasible solution to (5) and let Aj, g=1,2,...,r be 
a partition of the interval Ag = (0,1). Also, let AZ C A be a set of arcs (i, 7) 
such that for all A € Ag, it satisfies 7;(A) — 7; (A) = ci; (A), that is, 
Az = {(i,9) © A: mi(A) — my(A) = cig (A), AE Ag}. 
For the arcs not in the set AZ, assuming A € Ag, if m — 1; < cj;(A), then 
Lig = 0 and Vij = 0, and if Ti 1 > Cij(A), then Lig = Wij and Vig = 1 or 
xij =0 and y;; = 0. We define sets Lg and U, for q=1,2,...,r, as follows: 


Lg={,9) € Asm — 17 < cy(A), A € Ag}, 


Ue={G.7) € Ate ay > GglAly AEA, 


For each interval A,, we now define the restricted primal problem as 
follows. This problem tries to find a feasible solution for the problem of the 
parametric single-objective minimum cost-time network flow (4) using the the 
arcs in the set AZ, only. Consider the following restricted primal (RP(A7 ))) 
problem: 


40 O. Baghani and 8. Ghafoori 


(RP(A;)) min z = » Sk, 


keEV\{1} 
s.t. S- 15 S- Via — $y + So (-1)* sx = bi, (7) 
(L,j)EAZ (i,1)E AZ k=2 
S- Lkj — ain — (—1)** 8x = bg Vk € V\{I}, 
(k,j)EAG (i,k)EAP 
Lig S Wag hag VG,5) € Ay, 
big 20, vig € {0,1} V(i,5) € Ay, 
8 20 VREEY, 
where ; 
r™ Lo by <0, 


in which for each k € Y, by is defined as 


by, = by + Me Lik — S- Chi: 


(i,k)EU, (k,i)EUg 


Problem (7) has a feasible solution as follows: 
vij=0, yig=0, 81=0, se=lbx| (k © V\{1}). 


Note that RP(A7) does not depend on A. Assume that (2, 9, §) is the optimal 
solution of problem (7). If the optimal value of the problem RP(A7), 2, is 
zero, then the optimal solutions for problem (4) are as aj; = %ij and yj; = Ji; 
for all (7,7) € Aj, and xj; = 0 and yj; = 0 or aj; = uy and yj; = 1 for 
remaining arcs in A. Now if Z > 0, then the given solution is not a feasible 
solution to the primal problem (4). At this time, either a new feasible dual 
solution should be found that improves the objective function (7) or it can be 
concluded that the primal problem is infeasible. Since (#, 9, 8) is an optimal 
solution to problem (7) with the arcs in A, so a new arc in A\AZ is needed 
to reduce the value of the objective function in problem (7). We now consider 
the dual of the RP(.A7) as follows: 
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(DRP(AZ)) max So bete—- So Sis, 
kev (i, j)EAZ 
s.t. Ti — yj — big <0 for all (i,7) € Aj, 
[ij tig — O57 <0 for all (7,7) € Ags 
Tm >O0 (8) 
(—1)**r, — (—1)** x, <1 for all k € V\{1}, 
i050 for all (i,j) € Az. 


Suppose that (7, fi, 6) is an optimal solution for DRP(.A7 ). According to the 
complementary conditions, we obtain 


i) (iA; — sg) Biy = 0 for all i,j € V & for all (i,j) € Az, 


i) (- fig big — 6.) big = =0 for all (7,7) € Az: 
iti) (-4 y+ wig ig ) ig = =0 for all (i,7) € Aj, 
iv) (i- is) iy = 0 for all (i,j) € A. 


If0< Lij < Uay; then [hig = 0 and 7; — Tj = 0. In addition, if 7; — Tj <0, 
then Li; = 0 and Vig = 0. If 7; +9; > 0, then fig > 0 and Li; = Dig Wij - From 


the dual constraint fijjuij — bi < 0, it can be deduced that big > 0, therefore 


Vij = 1 and so Lij = 0. 
We now define the following sets: 
Ay” ={i,9) © AZ: ti — 7; = OF, 
Ly = {(i,5) € AZ : ti — tj < 0} = {Gi 7) € AT\ AT 
Ur ={(i,5) € Ag : ti — tj > OF} = {(t,9) € AT \AT 


1}. 


Q 


+ Lig = Vig, Yig = 


Now Suppose that (%, fi,5) is a new solution to the dual problem (5), 
which is defined as follows: 


nA) =a +O(A)z, 
f(A) = w+ OA)A, (9) 
(A) = 5 + A(A)4, 


where (7, 4,6) is a dual feasible solution for the original dual problem (5), 
(7, fi, 6) is a dual feasible solution for DRP(A7), g=1,.--,7r, and @(\) > 0 
For (i,7) € Aj, the new solution (7, /4, 5) is a dual feasible solution for (5), 
because 
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Ty = Tj — frig -— izi= (mi + O(r a) = (7 + (A); ) = (His + (A) jus) = Ci (A) 


= (n= m5 — Hy —6400)) +000) (m1 2s ~ Hu) 
< 0. 


For (i,j) ¢ AZ, the dual feasibility depends on the parameter (A). The 
new arc (i,j) € A\A7, which will be added to the current set AT, must 
reduce the value of objective function RP(A7). To this end, any are (i, 7) € 
A\ Aj satisfying either (7,7) € Lg (i-e., mi — mj < cij(A)) and 7 — 7; > 0 or 
(7,9) € Ug (ie., mi — mj > cig(A)) and 7; — 7%; < 0 can be included to AT 
provided that the dual feasibility for problem (5) is valid. Thus for keeping 
this condition, (A) is chosen as 


Cig (A) = Ti + HG + May 
Ty — 15 — frig 


((,5) € Uy, i — Ay <0)}. 


6(A) = min { (GA) € Ly, ti— i; > 0), 


(10) 


This process continues until 2 = 0. At this time, we will have an optimal 
solution for the primal problem (4). If the condition 2 = 0 is not satisfied, 
then an arc such as (i,j) € Lg such that 7; — 7; > 0 or (é,7) € Uy so that 
iti; —itj; <0 cannot be found. Therefore, the starting problem (4) is infeasible. 


The primal-dual algorithm for solving the BOMCTNFP is summarized 
in Algorithm 1. 


5 A biobjective network flow example 


Consider the following cities network. where 9 cities of Iran with 17 arcs 
are shown in Figure 1 (each city is a node). For each arc (7,7), there is a 
triple (c;;, tij, wiz) specifying the cost (million rials), time (hours) and ca- 
pacity (million tons), respectively. Also, for each city i, there is a scalar b; 
that indicates the value of supply or demand (million tons). As can be seen 
from Figure 1, three cities produce a special commodity, that is, Mashhad, 
Isfahan, and Ahwaz, and the rest cities are consumers. We trade 13.5 million 
tons of all network goods, and supply all demands with minimum cost at the 
minimum time. For simplicity, each city of the network is labeled with the 
numbers 1 to 9, as Table 1. 
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Algorithm 1 Primal-dual biobjective minimum cost-time network flow algorithm 


1: 


Set Ao := (0,1), and let (7°, u°, 6°)(A) = (0,0,0), A € Ao be an initial 
feasible solution of (5); 

Set C := {lp}, where lo := (Ao, 45 Sie. 3S 
0, (n°, n°, 8°)(A)); 

While C 40 do 


Select |, ¢ C and solve RP(A7) with 1,; 
(a) If the optimal value is 0, then any optimal solution is an 


efficient solution to (3). Set C := C\{lq}. 


(b) Else solve DRP(A7). Let (7, f,5)(A), A € Ag be an optimal 
solution for DRP(A7); 

i. If there is no arc (7,7) such that 7; — 71; < cj;(A) and 
i — Tj >0or 

MT, — 1; > Gj(A) and 7; — 7%; < 0, then (3) is infeasible 
and STOP. 

ii. Else set Aga), q’ =1...,7’ of Ag employing (10) and 
compute the new feasible solution (1 , 27 , 597 )(A), 
q =1,...,7r’ according to (10). 

Put C :=C\{Iq} Utlaq’}- 
End while 


[ oa %, 
Tabriz ~ 
4 TenTaN ae: Mashhad 
ey (22) 


G@) fcuth,(r) 


Zahedan 4, 
(100,9,1) %, 


1), sada abbas 


Figure 1: Graph of 9 cities of Iran with 17 arcs 


Keshavarz and Toloo [16] computed the supported efficient solutions and 


supported nondominated points for only two values of the scalarization pa- 
rameter A = 0.05, 0.5, which is reported in Table 2. Utilizing the primal-dual 
technique proposed in this paper, we have obtained all the supported efficient 
solutions and nondominated points of the problem that are introduced in Ta- 
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Table 1: Labelling the cities of the network 


City | Tabriz Tehran Esfahan Ahvaz Yazd Shiraz Mashhad Zahedan Bandar abbas 
label i 2 3 4 5 6 i 8 9 


Table 2: Efficient solutions and nondominated points proposed in [16] 


» Supported efficient solution Supported nondominated point 
X39 = 3.2, 236 = 2.1, 23, = 1.8 
0.05 a7, = 1.4, 2%. = 0.8,275, = 0.4 (cx, ty) = (1270.9, 66) 


Lag = 2.2,25, = 2.0, egg = 2.1 


i, =2804. S014, =12 


0.50 23, = 3.0, 275 = 1.4, 272 = 0.8 (cx, ty) = (1242.5, 77) 
Lig = 3.0, 27, = 1.2, 2G = 0.9 
Leg = 1.2 


ble 3. As can be seen, some new supported efficient solutions are computed 
that such solutions were not investigated earlier in [16]. These new efficient 
solutions supply a decision maker with more alternatives, who prefers less to 
more in each objective. 
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Table 3: Efficient solutions and nondominated points computed by our method for the 
network of Figure 1 


Supported efficient solutions 


Lig Vi,j Supported nondominated points 
U5 = 3.2, 03, = 2.1,73, = 1.8 ¥30 = 136 = 1,935 = 1 
Lrg = 1.4, 272 = 0.8, 25, = 0.4 ro = 1, y7g = 1,43, =1 
Lig = 2.2,25, = 2.0, 2% = 2.1 Yag = 14, = 1, ¥§ = 1 (cx, ty) = (1270.9, 66) 
15 = 13, = C43 = T3_ = 0 Y25 = ¥31 = Ya3 = ¥56 = 0 
Lig = 59 = C75 = UE9 = 0 Y¥5g = Y50 = Y75 = Yao = 9 
U3q = 2.8, 23, = 0.1, 23, = 1.2 Yzo = 1,936 = 1. y3, =1 
£35 = 3.0, 275 = 1.4, 27, = 0.8 ¥35 = 1y7g =1yzg =1 
Lig = 3.0,277, = 1.2, 2% = 0.9 Yig = 1 ya, = 1 vg =1 
tig = 1.2 Yeo = 1 (cx, ty) = (1242.5, 77) 
131 = 135 = €43 = T36 = 0 Y1 = Yo5 = Yaz = Y56 = 0 
Tig = C75 = CZo = 0 ¥53 = V75 = Ygo = 0 
£3, = 0.4, 235 = 2.8, 23, = 2.9 ¥31 = 1, y32 = 1,435 =1 
£36 = 3.0,23, = 2.0, 273 = 2.0 ¥36 = 1,y4, = Lag =1 
Lig = 0.2, erg = 1.1, 2G = 1.0 Yig = 1, yk = 1, vg =1 
Lrg = 1.4, 2%, = 0.8 Yro =1,y7g =1 (cx, ty) = (1381.4, 85) 
131 = T35 = U55 = 0 Yd1 = Yo5 = 56 = 0 
©5g = X75 = Tg9 = 0 ¥58 = Y75 = Yeo = 0 
£3, = 0.4,235 = 2.8,735 = 1.8 ¥31 = 1, y32 = 1,y35 =1 
3g = 21,23, = 2.0, 03, = 2.2 ¥36 = 1 yd, = Lyag =1 
Lég = 2.1, 279 = 1.4, 27 = 0.8 ¥éo = 1, y7g = 1yzg = 1 (cx, ty) = (1250.9, 70) 
131 = 135 = L43 = T3_ = 0 Y21 = Yo5 = Yaz = Y56 = 0 
@5g = 59 = C75 = Teo = 0 Y¥53 = Y50 = Y75 = Yso = 0 
£3, = 0.4, 235 = 2.8, 23, = 3.0 ¥31 = 1, y32 = 1,435 =1 
£36 — 0.9,23, = 2.0, 236 = 2.2 ¥36 = 1,44, = 1yuag =1 
Leg = 1.2, 2%) = 0.9, 27, = 1.4 Yeo = 1, yé9 = Lyre =1 (cx, ty) = (1156.9, 67) 
131 = 135 = Tj3 = T3_ = 0 Y21 = Yo5 = Yaz = ¥56 = 0 
©5g = 75 = Tig = Teo = 0 Y¥58 = Y75 = Vis = Yao = 0 
v3, = 1.2, 2735 = 2.8, 23, = 1.8 ¥31 = 1,y¥32 = 1y3, =1 
23g = 1.3, 23, = 1.2, 23, = 3.0 ¥36 = 1, y4, = 1yag =1 
Leg = 21,27, = 1.4, 27. = 0.8 Ye9 = LYfg = 1 yFg = 1 (cx, ty) = (1248.5, 70) 
131 = ©}5 = T43 = Tq = 0 Y21 = ¥o5 = Y43 = ¥56 = 0 
T5g = hq = C75 = Taq = 0 ¥5g = Y50 = Y75 = Yso = 0 


6 Computational experiments 


The primal-dual approach for solving five instances of the BOMCTNFP has 
been implemented using MATLAB R2017a software, on a laptop equipped 
with Intel Core i5-8250U 1.80GHz processor and 8 GB of RAM, to find all 
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their efficient extreme points. Table 4 shows that when the number of nodes 
and arcs grows as compared to the number of supported efficient solutions, 
the iterations and computational cost of the algorithm extremely increase. 


Table 4: Numerical results by Algorithm 1 for several instances of BOMCTNFP (nses 
stands for the number of supported efficient solutions) 


The parameters of each instance 


Instances Nodes Arcs Sources Sinks max cj max ;,; Wij nses_ Iterations CPU time (sec.) 


il 4 6 2 2 4 2 Li 1 2 0.02 
2 5 8 3 2 7 3 1-3 2 5 0.98 
3 8 14 3 5 26 3 1-3 3 43 1.24 
4 9 16 5 4 34 4 1-4 3 137 2.78 
5 14 29 10 4 62 14 1-5 6 198 4.53 


7 Conclusion 


The purpose of this paper was to solve the mixed biobjective problem of 
minimum cost-time network flow. To this end, we first formulated the mixed 
biobjective problem and proposed the primal-dual algorithm for solving it. In 
this article, an example with 9 city nodes and 17 arcs was solved. We found 
6 supported efficient solutions in the decision space. In comparison with 
the method presented in [16], we could find some new supported efficient 
solutions to provide a decision maker with more alternatives. 
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